Trapping a particle of a quantum walk on the line 



Antoni Wojcik, 1 Tomasz Luczak, 2 Pawel Kurzynski, 3 ' Andrzej 
Grudka, 1 Tomasz Gdala, 2 and Malgorzata Bednarska-Bzd§ga 2 

1 Faculty of Physics, Adam Mickiewicz University, Umultowska 85, 61-614 Poznari, Poland. 
2 Faculty of Mathematics and Computer Science, Adam Mickiewicz University, Umultowska 87, 61-614 Poznari, Poland. 
3 Centre for Quantum Technologies, National University of Singapore, 3 Science Drive 2, 117543 Singapore, Singapore 

(Dated: December 7, 2011) 

We observe that changing a phase at a single point in a discrete quantum walk results in a rather 
surprising localization effect. For certain values of this phase change the possibility of localization 
strongly depends on the internal coin-state of the walker. 
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I. INTRODUCTION 



Quantum analogs of random walks have been intro- 
duced in hope to device efficient algorithms of quan- 
tum computations [T] and to study new quantum effects 
[2J. For the last decade they have attracted much atten- 
tion and a number of non-classical properties of quantum 
walks have been described [3HS] and experimentally ver- 
ified [6 15 . In particular, an enhancement and quantum 
suppression typical for chaotic quantum systems have 
been predicted and observed for different variants of a 
one-dimensional quantum walk |16H18j . 

Here, we describe yet another feature of a quantum 
walk on an infinite chain: a sharp allocation of a particle 
around a distinguished node, caused by a local modifica- 
tion of the phase of a particle. We remark that allocations 
in quantum walks were observed before (see e.g. [T5lJ29| ). 
In particular, related research include studies of localiza- 
tion in coined quantum walks with one defect coin at 
origin [29-31 . Here we consider a simple inhomogeneity 
caused by a single phase defect at origin, which neverthe- 
less leads to interesting localization effects and, possibly, 
can be observed experimentally. We find that the num- 
ber of localized stationary states strongly depends on the 
phase defect. We give analytical formula for these states 
and show that for some phase defects the effect of local- 
ization can be controlled via coin states. 

Before we proceed, let us recall the standard model of 
a coined quantum walk. The evolution of the system is 
given by a unitary operator U defined on a tensor product 
of two Hilbert spaces H c <8> H n , where H c is spanned 
by the coin states |c), (c = 0, 1), while the basis of H n 
consists of vectors |n), (n E Z), which correspond to the 
position of a particle on the chain. In a single step of the 
evolution of the system we toss a coin, i.e., change the 
coin state |c), and then move the particle to one of the 
two neighboring states of the chain determined by |c). 



The corresponding unitary operator is given by 

U= C£2\c){c\®S c )(H®I), (1) 

c 

where 



is a Hadamard matrix and S c is a translation operator 
defined as S c \n) = \n + (— l) c+1 ). The state vector 

l*) = £(«»|o>l»>+/gi>N) (3) 

n 

of the system evolves in (discrete) time t according 
to |\£(t)) = f7*|*(0)). Because the evolution of the 
state supported on even positions does not interfere with 
the evolution of the state supported on odd positions, 
therefore it is convenient to consider double steps, i.e., 
|\P(t)) = Cf 2t |*(0)). In this case the walk that is initially 
localized in position n = occupies only even positions. 

II. SINGLE-POINT PHASE DEFECT AND 
LOCALIZATION 

Let us now break the translational symmetry of the 
above model by modifying the phase of the particle by 
lu = oj(4>) = e 2m< ^ [(f) £ (0, 1)) each time it passes through 
node n = 0. More precisely, we replace U by U$ = 
(EJ C K C I ® S c)(C ® I), where Sf\n) = e 27ri ^.°|n + 
( — l) c+1 ) and <5 rii o is Kronecker delta. We note that if the 
state of the coin is measured after each step, the model 
corresponds precisely to the classical random walk. 

It turns out that the above symmetry breaking at node 
zero results in an emergence of new localized eigenvectors 
of the corresponding quantum walk. These eigenvectors 
and the corresponding eigenvalues can be found by solv- 
ing a set of recurrence equations generated by operator 
U^. We have "atypical" equations — these ones which in- 
volve amplitudes of probability corresponding to particle 
at node 



'Corresponding author. E-mail: cqtpkk@nus.edu. sg 



\/2a_i(t + l) = Luao(t) + ujf3o(t), 
y/20x(t + l)=wa o (t)-u0o(t) (4) 



2 



and "typical" equations, i.e., the ones whose right hand 
side does not refer to amplitudes corresponding to parti- 
cle at node 



V2a n (t + 1) = a n+ i(t) + P n+ x{t), 
V2Pn{t + 1) = a„-i(t) - p n -i(t). 



(5) 



In the above we denoted component of a state vector 
corresponding to particle in node n and coin in state 
|0) at time t by a n (t) and component corresponding to 
node n and coin in state |1) by (3 n (t). For the double step 
operator [/J and for stationary states these equations are 
time independent and take the form 

2Aao = <5i + /?i + loo>o — ujpo, 
2X(3 = ujao + w/3 — a_i + P-i, 
2Aa_x = ua + u)j3 a + — f3_i, 

2\(3i = ai + Pi - wd + o;/3 (6) 

and 

2Xa n = a n+ i + p n+ i +a n - p n for n ^ 0, -1, 
2\j3 n = a n + p n - a„_i + j3 n -i for n ^ 0, 1, (7) 

where A denotes an eigenvalue of C/J and for convenience 
we substituted a n = oc2n and /3 n = fan- 

The solution to Eqs. [6] and [7] (for detailed derivation 
see the Appendix) yields 



A+ = 



u) — 2uj 2 + uj 3 ± iuj(l — uj + oj 2 ) 
1 - 2u + 2uj 2 : 



(8) 



whereas the components of the localized stationary state 
vectors which do not correspond to the particle at node 
are given by 

=Cx r \ 



and 



where 



i for n > 1, 

= C(u> qp «w ± for n < — 1 

= C(l — uj =F iu;)x5 for n > 1, 

= ^iCx± n for n < -1, 



x± = 



2 cos(2tt</>) =f 2 sin(27T0) - 3 ' 



(9) 



(10) 



(11) 



The components of the localized stationary state vectors 
which correspond to the particle at node are 



* (±) - 
'o — 

?(±) _ 



C. 



In the above C is a normalization constant 



C = 



l+x± 



(12) 



(13) 



It follows that there are at most two stationary states to 
which we refer as to \Q±). 
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FIG. 1: The plots of absolute values of x+ (dashed) and 
x- (solid) as a function of <j>. The shaded region denotes the 



allowed values of x± . 
for<^G (0,|). 



X- 



< 1 for 4> £ (t, l) whereas |a;+| < 1 



III. DISCUSSION 



First, let us observe that the number of stationary 
states depends on 4>. Due to the normalization condi- 
tion both a n and f3 n should tend to zero as n tends to 
±oo. Therefore, the parameter x± defined in (111 has 
to obey |a;±J< 1. This is true only for certain values of 
</> (see Fig. Ill), namely < 1 for <fi 6 (|, l) whereas 
|x_|_| < 1 for <fi € (0, 1). Therefore, in the range (j> £ (0 
there is only one localized stationary state |$ + ) 
range </> g ( j, f ) there are two such states |$ + ) and 1$-), 
and in the range <fi € [|, 1) there is again one localized 
stationary state, but this time it is \$-)- 



in the 



Next, let us study the localization properties of the 
quantum walk with single-point phase defect. We say 
that the walk exhibits localization if in the limit of infi- 
nite (even) time the probability of occupying the initial 
position does not tend to zero |30j . From the previous 
paragraph we know that there are at most two normal- 
izable eigenstates of Z7| that correspond to stationary 
states |$±). By the analogy to textbook problems of 
a quantum particle in the presence of a potential trap, 
these states can be thought of as bound states. The re- 
maining states correspond to scattering states. We ex- 
pect that an initially localized state whose decomposition 
into eigenstates of the evolution operator does not con- 
tain bound states does not exhibit localization (see [30] 
and references therein) . Therefore, the overlap of the ini- 
tial state of the walk with bound eigenstates can be used 
to measure the localization properties of the walk. Note 
that this overlap, which we denote as F, is definitely a 
lower bound for the probability of occupation of the ini- 
tial position in the limit of infinite time. 



For initial state of the form \ipo) = a\< 
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FIG. 2: The plots of overlap F of initial states of the quan- 
tum walk and the localized stationary states of 1/2. The 
solid line corresponds to the initial coin state |0), whereas the 
dashed line corresponds to the initial coin state 775 (|0) 



/3|l)|n = 0), the total overlap is: 

F = F + for G (0, \] 



F = F., 

F = F- for 



F_ for 6 (|, |; 
[f,l). 



where 



F± = K*±|^o>| 2 = 



l + a;± 



(14) 



(15) 



Fig. [2] shows the plots of F for initial coin states |0) 
and TygdO) + i|l))< This example evidently shows that 
localization properties depend on the initial coin state of 
the particle. In particular, there are two initial states 
= 4g(|0)|n = 0) ± i\l)\n = 0)) that are orthogo- 
nal to localized stationary states |$^), respectively. As a 
result, we observe that the initial state \il>o) exhibits lo- 
calization only in the range <f> € (0, |), whereas the state 
\iPq) exhibits localization only in the range (f> 6 (|>l)- 
In the Fig. [3] we present probability distributions for a 
quantum walk with phase defect at origin corresponding 
to (f> = g for initial states \iPq) for which both localization 
and the lack of localization are well visible. 



IV. EXPERIMENTAL PROPOSAL 

Here we propose an experiment to verify our theoreti- 
cal results. It is a simple modification of the experiment 
performed by Broome et. al. [T3] that employs a single 
photon "walking" in an array of beam splitters (beam 
displacers) where the polarization of the photon plays 
the role of the coin. We believe that our proposal can 
be relatively easy accommodated to any other quantum 
walk implementation. A pair of photons is created in a 
process of spontaneous parametric down-conversion. The 
presence of the photon that is going to be injected into an 



FIG. 3: Probability distribution after 50 steps of a quantum 
walk with a single phase defect at origin corresponding to 
<j> = |. Solid line — lack of localization for the initial coin 
state 7^75 ( 1 0) — Dashed line — localization at origin for 

the initial coin state 7^75 ( 1 0} + i|l)). 
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FIG. 4: a) Experimental proposal for four steps of a quan- 
tum walk with a phase defect at position n = 0. It em- 
ploys a single photon in an array of optical elements and is 
based on a setup introduced in |13| . Bold (blue) extended 
rectangles represent calcite beam displacers which realize a 
coin-dependent position shift, thin (red) extended rectangles 
represent half- wave plates/quarter- wave plates which realize 
coin operations, and small (yellow) rectangles represent phase 
shifters which implement a phase shift at position n — 0. Be- 
low we present expected probability distributions after four 
steps for a phase shift corresponding to <f> = | . b) Initial coin 
state 775OO) c) Initial coin state tt=(|0) 



array of calcite beam displacers and wave plates is her- 
alded by a detection event of its twin. Before the walk 
starts, the initial state of the polarization (coin) is pre- 
pared via application of a proper wave plate. The scheme 
of an array implementing a quantum walk with a single 
phase defect at position n — is presented in Fig. [4] 
Calcite crystals displace the path of the photon with re- 



4 



spect to its polarization state, therefore they are capable 
of realizing a coin-dependent position shift. Wave plates 
change the polarization state, hence they realize the coin 
operation. Finally, properly placed phase shifters can 
implement a phase shift at origin. 

Due to the fact that in realistic experiments one always 
encounters inevitable errors it is of great importance to 
determine the conditions under which the investigated 
phenomena is observable. It is obvious that the more 
steps of a quantum walk one wants to implement, the 
more robust to noise the experimental setup has to be. 
Fortunately for us, the first signs of the coin-dependent 
localization can be observed after only four steps of the 
corresponding walk (see Fig. |4]b and c). Note that the 
authors of [13] implemented six steps. 

An interesting addition to the experiment would be to 
investigate the effects of decoherence on localization. It is 
well known that decoherence transforms quantum walks 
into classical random walks |32] for which interference ef- 
fects play no role, hence the phase shift at position n = 
would cause no effect and the corresponding walk should 
exhibit Gaussian pattern irrespective of the initial coin 
state. Note that tunable decehorence is within experi- 
mental reach since the authors of [13] implemented it via 
change of relative angles between calcite crystals. 



V. FURTHER RESEARCH 

The state-dependence of localization can be used to fil- 
ter and to trap particles of discrete time quantum walks 
with respect to their internal coin state. This effect can 
find applications in quantum algorithms and in quantum 
state engineering. Another interesting perspective is to 
investigate localization due to a single phase defect in 
quantum walks in two (and more) dimensions. In partic- 
ular, it would be interesting to investigate single-point lo- 
calization on grids and honeycomb lattices, since the for- 
mer represent natural generalization of a one-dimensional 
chain and the latter are interesting from the point of view 
of the recent intensive research on the material properties 
of graphene. 
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A. Appendix 

Here we give detailed solution to Eqs. Q and Q. 
Recall that we use the following substitution a n = a^ n 
and f} n — /3 2n - From Eqs. Q one obtains 

(3 n +i = ("n+i - Aa w ) for n ^ -1, (16) 

Substituting the above into the second equation in ^ 
we have 

Xa n+ i — 2(A 2 — A + l)a„ + Aa„_i = 

for n 0,1 (17) 

The general solution to this recurrence equation is of the 
form 

a n = C+x n + C-x~ n , 



5 



where x as well as i 1 satisfy equation 



Xx 2 



2(A 2 - A + l)x + A = for i ^ -1, 0, 1 (18) 



and C + , C_ are constants. Due to normalization con- 
straint for n — > ±oo we must have \a n \ — ► 0, therefore 



a n — C + x n for n > 1, 
i„ = C_a;~ n for n < —1, 



(19) 



where x is the solution to (18) that satisfies |x| < 1 



Substituting Eqs. (j 1 9|) into Eq. (j 1 6|) we have 



a 



Pn = C_ 



x— A ™n— 1 
+ A-l X 

l-Az -t 
A-l x 



for n > 2, 
for n < — 1. 



(20) 



From the first and the fourth equation of ^ we obtain 



a 



* -A)). 



(21) 



Substituting Eqs. (19) and (20) into the above we have 
a = C+. (22) 



Observe that due to the above the first equation in ( 20 ) 
is also valid for n = 1. From the second and the third 
equation of Q we obtain 



A) = ^09-i+a-i(A-l)). 



From Eqs. (19), (20 1 and (18) it follows that 



% = C. 



1 - Air 
A - 1 ' 



(23) 



(24) 



From Eqs. ( |T9[ ) and p0[ ) and from the first and the sec- 
ond equation of (JsJ) we have 

ujp = C + (x + uj-2\+fE 2 [), 
C + u> = P (2\-x-u + x^), (25) 

which together give 



ui 2 = (x + oj — 2A - 



x A . , A 1. . . . 

-)(2A -X-W + x- — ).(26) 



A-l 



1-Xx' 



Taking into account Eq. ( 18 ) the solution to the above 
equation yields 



oj - 2u> 2 + lo 3 ± -UJ + LO 2 ) 

1 - 2uj + 2uj 2 



and 



x± 



lo 



uj 2 -3u + l±i(uj 2 -1)' 
which can be expressed as 

1 

X± ~ 2cos(27rai>) =F 2sin(27T(A)^3' 
From Eqs. pi}, p5|, ( 27 1 and (29D one 



(27) 



(28) 



one finds 



C_ = C + (uj =F ± i), 



(29) 



(30) 



where the ± sign depends on the signs in Eqs. (27) and 
(29). As a consequence, we arrive at 



fl(±) 

rn 



Cxi for n>l, 

C(lu =f iuj ± i)x± n for n < — 1, 
= C(l — w =F «w)a;" for n > 1, 



and 



= T*Cx7 n 



for n < — 1 



(31) 



(±) 



(32) 



where we substituted C = C+. The normalization con- 
stant C can be easily evaluated since the sum of proba- 
bilities \a n \ 2 and \(3 n \ 2 gives rise to the sum of geometric 
series. It yields 



C 



1 + £C± 



(33) 



